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Abstract 

Elaborating on M. Pavon, J. Math. Phys. 40 (1999), 5565-5577, we de- 
velop a simplified version of a variational principle within Nelson stochas- 
tic mechanics that produces the von Neumann wave packet reduction 
after a position measurement. This stochastic control problem parallels, 
with a different kinematics, the problem of the Schrodinger bridge. This 
gives a profound meaning to what was observed by Schrodinger in 1931 
concerning Schrodinger bridges: "Merkwiirdige Analogien zur Quanten- 
mechanik, die mir sehr des Hindenkens wert erscheinen" . 



1 Introduction: Schrodinger's problem 

In 1931/32 P3I2], Schrodinger considered the following problem. A cloud 
of N Brownian particles in M n has been observed having at time to an 
empirical distribution approximately equal to po(x)dx. At some later 
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time ti, an empirical distribution approximately equal to pi(x)dx is ob- 
served. Suppose that p\{x) considerably differs from what it should be 
according to the law of large numbers (N is large), namely 

rti 

/ P(t ,y,t 1 ,x)p (y)dy, 

where 

n 

p(s,y,t,x) = [2ir(t - s)) 2 exp 

is the transition density of the Wiener process. It is apparent that the 
particles have been transported in an unlikely way. But of the many 
unlikely ways in which this could have happened, which one is the most 
likely? 

In modern probabilistic language, this is a problem of large devia- 
tions of the empirical distribution PJ. By discretization and passage to 
the limit, Schrodinger computed the most likely intermediate empirical 
distribution as iV — > oo. It turned out that the optimal random evolu- 
tion, the Schrodinger bridge from po to p\ over Brownian motion, had at 
each time a density p(-,t) that factored as p(x,t) = <p(x,t)<p(x,t), where 
<p and <fi are a p-harmonic and a p-coharmonic functions, respectively. 
That is 

0(t, x ) = J pit, x, t u y)(p(t u y)dy, (1) 
4>(t, x) = J pfo, y, t, x)(j)(t , y)dy. (2) 

The existence and uniqueness of a pair ((/>,</>) satisfying (HJ)-© and the 
boundary conditions (J)(x,to)(fi(x,to) = po(x), (f>(x,ti)(f>(x,ti) = pi(x) was 
guessed by Schrodinger on the basis of his intuition. He was later shown 
to be quite right in various degrees of generality by Fortet [4 , Beurlin 
Jamison jB], Follmer 0. Jamison showed, in particular, that the 
Schrodinger bridge is the unique Markov process {x(t)} in the class of 
reciprocal processes (one-dimensional Markov fields) introduced by Bern- 
stein |7j having as two-sided transition density 

p(s,x;t,y)p(t,y;u,z) 

q(s,x;t,y;u,z) = r , s < t < u, 

p(s, x; u, z) 



\x — y\ 
2(t-s) 



s <t, 
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namely q(s,x;t,y;u, z)dy is the probability of finding the process x in 
the volume dy at time t given that x(s) = x and x{u) = z. Schrodinger 
was struck by the following remarkable property of the solution: The 
Schrodinger bridge from p% to po over Brownian motion is just the time 
reversal of the Schrodinger bridge from po to p%. In Schrodinger's words: 
"Abnormal states have arisen with high probability by an exact time 
reversal of a proper diffusion process" . This led him to entitle : "On 
the reversal of natural laws" A few years later, Kolmogorov wrote a paper 
on the subject with a very similar title Moreover, the fact that the 
Schrodinger bridge has density p(x,t) = (f)(x,t)(f)(x,t) resembles the fact 
that in quantum mechanics the density may be expressed as p(x, t) = 
ijj(x,t)^)(x,t). This analogy has inspired various attempts to construct 
a stochastic reformulation of quantum mechanics starting from 

[TJ 12 El ■ In order to discuss a more general Schrodinger bridge problem, 
we recall in the next session some essential facts on the kinematics of 
finite-energy diffusions as presented in [T3J EJ E3 EE] . 

2 Elements of Nelson-Follmer kinematics 
of finite energy diffusions 

Let (fl, T, P) be a complete probability space. A stochastic process 
< t < ti} mapping [to,^i] into L^(Q, J 7 , P) is called a finite- 
energy diffusion with constant diffusion coefficient I n a 2 if the path £(ct>) 
belongs a.s. to C([t Q , ti}; R n ) (n-dimensional continuous functions) and 

£(*)-£(*)= /* P(r)dr + aw+(s,t), t <s<t<t u (3) 

J s 

where the forward drift f3{t) is at each time t a measurable function of 
the past {£(r);0 < r < t}, and is a standard, n-dimensional 

Wiener difference process with the property that w + (s, t) is independent 
of {£(t); < r < s}. Moreover, f3 must satisfy the finite-energy condition 

E j^* 1 0{t) ■ /9(r)dr} < oo. (4) 

We recall the characterizing properties of the n-dimensional Wiener dif- 
ference process w + (s,t), see O Chapter 11] and [TSJ Section 1]. It is a 
process such that w+(t,s) = —w + (s,t), w + (s,u) + w+(u,t) = w + (s,t), 
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and that w + (s,t) is Gaussian distributed with mean zero and variance 
I n \s — t\. Moreover, (the components of) w + (s,t) and w + (u,v) are in- 
dependent whenever [s,t] and [u, v] don't overlap. Of course, w + (t) := 
w + (t , t) is a standard Wiener process such that w + (s, t) = w + (t)—w + (s). 
In [T3], Follmer has shown that a finite-energy diffusion also admits a 
reverse-time differential. Namely, there exists a measurable function *~y(t) 
of the future {£(t); t < t < ti} called backward drift, and another Wiener 
difference process u>_ such that 

-£(s) = f\(r)dT + aw^(s,t), t < s < t < t x . (5) 
Moreover, 7 satisfies 

rti 



E 



jjf 7(r) -7(r)dT| < 00, (6) 

and is independent of {£(V);i < r < tx}. Let us agree that dt 

always indicates a strictly positive variable. For any function / defined 
on [i , ti], let 

d+ffl = f(t + dt) - f(t) 
be the forward increment at time t, and 

d„f{t) = f{t)-f{t-dt) 

be the backward increment at time £. For a finite-energy diffusion, Follmer 
has also shown in ^3] that the forward and backward drifts may be ob- 
tained as Nelson's conditional derivatives, namely 



and 



the limits being taken in L^fl, J 7 , P). It was finally shown in that 
the one-time probability density p(-,t) of £(£) (which exists for every 
t > to) is absolutely continuous on M n and the following duality relation 
holds W > 

i?{/3(t)-7(t)K(t)} = a 2 Vlogp(e(t),t), a.s. (7) 
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Remark 1 It should be observed that in the study of reverse-time dif- 
ferentials of diffusion processes, initiated by Nelson in ^7] and Nagasawa 
in JHJj see ^Hl |2D] and references therein, important results have ob- 
tained by A. Linquist and G. Picci in the Gaussian case [2U 1221 without 
assumptions on the reverse-time differential. In particular, their results 
on Gauss-Markov processes have been crucial in order to develop a strong 
form of stochastic realization theory |2I]-j21] together with a variety of 
applications [21]- [30 . 

Corresponding to Q and (j3J) are two change of variables formulas. Let 
/ : M. n x [0,T] — > R be twice continuously differentiate with respect 
to the spatial variable and once with respect to time. Then, if £ is a 
finite-energy diffusion satisfying (j3J) and (J3J), we have 

/(£(*).*) - f {£(*),*) = + 000 • v + y A ) KfrMfr 

+ /Vv/(e(r),r).d + «; + (r), (8) 

«/ s 

/(£(*),*) - Ms), s) = f g (J- + 7 (r) -V-^a) m(r),r)dr 

+ j ( < ^V/(£(r),r)-d^_(r). (9) 

The stochastic integrals appearing in (jHJ) and Q are a (forward) Ito 
integral and a backward Ito integral, respectively, see [T5| for the details. 

3 Schrodinger bridges 

The solution to the Schrodinger problem can be obtained by solving a 
stochastic control problem. The Kullback-Leibler pseudo-distance be- 
tween two probability densities p(-) and q(-) is defined by 

H( P , q ):=J R Jo g P ^p(x)d X . 

This concept can be considerably generalized. Let Q := C([t , ti), M. n ) 
denote the family of n-dimensional continuous functions, let W x denote 
Wiener measure on Q starting at x, and let 

W:= f W x dx 



5 



be stationary Wiener measure. Let D be the family of distributions on 
Q that are equivalent to W. For Q, P e O, we define the relative entropy 
H(Q, P) of Q with respect to P as 

H(Q,P) = E Q [\og^]. 

It then follows from Girsanov's theorem that fOJ EH HI] 

H(Q,P) = H(q(t ),p(t )) + E Q 

= H(q(tx),p(ti)) + £ Q [ f - 7 P (t)] • [ 7 Q W - 7 F (*)H -(10) 



rtl lh Q (t)-i p (t)]-h Q (t)-7 p (t)]dt 

t 2. j 

Here q(to) is the marginal density of Q at to> and 7^ are the forward 
and the backward drifts of Q, respectively. Now let p an d Pi be two 
everywhere positive probability densities. Let D(p ,Pi) denote the set of 
distributions in D having the prescribed marginal densities at to and t\. 
Given ? G D, we consider the following problem: 

Minimize H(Q,P) over D(p ,pi). 

In view of (jlUJl . this is a stochastic control problem. It is connected 
through Sanov's theorem [3J HI2] to a problem of large deviations of 
the empirical distribution, according to Schrodinger original motivation. 
Namely, if X 1 , X 2 , ... is an i.i.d. sequence of random elements on Q with 
distribution P, then the sequence P n [- X^=i ^x* £ '] satisfies a large de- 
viation principle with rate function if (•, P). 

If there is at least one Q in D(p , p x ) such that H(Q, P) < 00, it may 
be shown that there exists a unique minimizer Q* in D(p ,pi) called 
the Schrodinger bridge from p to p\ over P. If (the coordinate process 
under) P is Markovian with forward drift field b+(x,t) and transition 
density p(a, x, r, y), then Q* is also Markovian with forward drift field 

bf (x, t) = b%{x, t)+V log <f>(x, t), 

where the (everywhere positive) function <fi solves together with another 
function (j) the system (JIJ)-(j2J) with boundary conditions 

(j)(x, t )4>(x, t ) = po(x), <f>(x, t-y)4>{x, ti) = pi(x). 
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Moreover, p(x,t) = (f)(x,t)(f>(x,t),Vt G [to>£i]- This result has been suit- 
ably extended to the case where P is non-Markovian in [31] . For a survey 
of the theory of Schrodinger bridges with an extended bibliography see 

Consider now the following simpler problem: We have a reference 
stochastic model P G D. We think of P as modeling the macroscopic 
evolution of a thermodynamic system. Suppose we observe at time t\ 
the (everywhere positive) density p\ different from the marginal density 
of P. Thus we need to solve the following optimization problem 

Minimize H(Q,P) over Q G D(pi). 

where O(pi) denotes the set of distributions in D having density pi at 
t\. Let us assume that H(pi,p(t\)) < oo. In view of (fTUj). this stochastic 
control problem can be trivially solved. The unique solution is given by 
the distribution Q* having backward drift j p {t) and marginal density 
Pi at time t\. Thus, the result of measurement at time t\ leads to the 
replacement of the stochastic model P with Q*. Notice that the backward 
drift 7 P (t) is perfectly preserved by this procedure. Symmetrically, if we 
were to change the initial distribution at time to, the procedure would 
preserve the forward drift (3 p {t). 

4 Elements of Nelson's stochastic mechan- 
ics 

Nelson's stochastic mechanics is a quantization procedure for classical 
dynamical systems based on diffusion processes. Following some early 
work by Feynes [33] and others, Nelson and Guerra elaborated a clean 
formulation starting from 1966 [SU EH IHI3 ; showing that the Schrodinger 
equation could be derived from a continuity type equation plus a Newton 
type law, provided one accepted a certain definition for the stochastic 
acceleration. In analogy to classical mechanics, the Newton-Nelson law 
was later shown to follow from a Hamilton-like stochastic variational 
principle [56*] E?T] . Other versions of the variational principle have been 
proposed in gHl 123 H3 El • 

Consider the case of a nonrelativistic particle of mass m. Let {ij)(x, t);t < 
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t < ti} be the solution of the Schrddinger equation 

dip ih 



dt 



'.ID 



such that 



||V^6LU*o,+oo). 



(12) 

This is Carlen's finite action condition. Under these hypotheses, the 
Nelson measure P G D may be constructed on path space, [I2],|III], (HE 
Chapter IV], and references therein. Namely, letting Q := C([to, ti], M n ) 
the n-dimensional continuous functions on [tQ,£i], under the probability 
measure P, the canonical coordinate process x(t,u) = u{t) is an n- 
dimensional Markov diffusion process {x(t);to <t< ti}, called Nelson's 
process, having (forward) Ito differential 



dx{t) 



— V(a^ogV>(x(*),t) + StogV>(x(t),t)) 
m 



dt+ \ —dw(t), (13) 
V m 



where w is a standard, n-dimensional Wiener process. Moreover, the 
probability density p(-,t) of x(t) satisfies 



p(x,t) 



{x,t)\ 2 , wte[t ,h]- 



(14) 



Following Nelson [121 EHj, for a finite-energy diffusion with stochastic 
differentials we define the current and osmotic drifts, respectively: 



v(t) 



u{t) 



(3{t) - 7 (t) 



Clearly v is similar to the classical velocity, whereas u is the velocity due 
to the "noise" which tends to zero when the diffusion coefficient a 2 tends 
to zero. In order to obtain a unique time-reversal invariant differential 
[4T)j . we take a complex linear combination of obtaining 



x(t) — x(s) 
a 



1 



- 7 (r) 



dr 



[(1 



2 ' v ' 2 
i)(w+(t) - w + (s)) + (1 + i)(w4t) - w4s))} 



Let us define the quantum drift 

1 -i 



v q (t) 



P(t) + ^±~f{t)=v(t)-iu(t), 



S 



and the quantum noise 



w q {t) 



1 



— i 



w+(t) + 



l + i 



2 



2 



Hence, 



x(t) — x(s) 



L 



v g (r)dT + a[w q (t) - w q (s)}. 



(15) 



This representation enjoys the time reversal invariance property. It has 
been crucial in order to develop a Lagrangian and a Hamiltonian dy- 
namics formalism in the context of Nelson's stochastic mechanics in 
1101 EH EE! • Notice that replacing ©-© with (fTHjl. we replace the pair 
of real drifts (v,u) by the unique complex-valued drift v — iu that tends 
correctly to v when the diffusion coefficient tends to zero. 

5 Quantum Schrodinger bridges 

We now consider the same problem as at the end of Section El We 
have a reference stochastic model P G D given by the Nelson measure 
on path space that has been constructed through a variational principle 
[SZllSHlEni- This Nelson process x = {x(t); to < t < t{\ has an associated 
solution {if)(x,t) to < t < t\} of the Schrodinger equation in the sense 
that the quantum drift of x is v g (t) = -^V log ■?/>(£(£), t) and the one-time 
density of x satisfies p(x, t) = \if)(x, t\ 2 . Suppose a position measurement 
at time t\ yields the probability density pi(x) ^ \vf)(x,ti)\ 2 . We need 
a suitable variational mechanism that, starting from (P,p\), produces 
the new stochastic model in D(pi). It is apparent that the variational 
problem of Section |H] is not suitable as it preserves the backward drift. 
Since in stochastic mechanics both differentials must be granted the same 
status, we need to change both drifts as little as possible given the new 
density pi at time t\. Thus, we employ the differential (JT3J), and consider 
the variational problem: 

Extremize on (x, v q ) G (B(pi) x V) the functional 




It • Mx{t),t)-v q {t))dt 



} 
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subject to: x has quantum drift (velocity) v q . 

Here v q (x,t) = -A- V log ip(x, t) is quantum drift of Nelson reference pro- 
cess, and D(pi) is family of finite-energy, M n -valued diffusions on [t ,ti] 
with diffusion coefficient — , and having marginal p\ at time t\. Moreover, 
V denotes the family of finite-energy, C n - valued stochastic processes 
on [toj^i]- Following the same variational analysis as in we get a 
Hamilton- Jacobi-Bellman type equation 

+ v q (x, t) ■ V<p(x, t) - 7^A^(x, t) = ^V(^(x, t) ■ Vif(x, t), (16) 

with terminal condition (p(x,ti) = | log j^tl) ' Then x £ O(pi) with 
quantum drift 

h 

v q (x(t),t) + — Vp(x(t),t) 
mi 

solves the extremization problem. Write ip(x, t\) = p(x, ti) a exp[lS(x, ti)], 
and let {ip(x,t}} be solution of Schrodinger equation (fTTj) on [t , ^i] with 
terminal condition 

~ i % 

ijj{x,ti) = pi(x)2 exp[-S(x,t{)]. 



Next, notice that for t £ [toi^i] 



_ + „ g(M) . V -_A 



-0/ V V p(iC, ti 



where v q (x,t) = \ogtp(x(t). It follows that (p(x,t) := log^(x,t) 
solves (fTB)l. and the corresponding quantum drift is 

v q (x(t),t) + —yip{x{t),t) = — V log 

Thus, new process after measurement at time t\ (quantum Schrodinger 
bridge) is just the Nelson process associated to another solution i/j of the 
same Schrodinger equation. Invariance of phase at t\ follows from the 
variational principle. 
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6 Collapse of the wavefunction 

Consider the case where measurement at time t\ only gives the informa- 
tion that x lies in subset D of configuration space M n of the system. The 
density pi(x) just after measurement is 



where p(x,t\) is density of Nelson reference process x at time t\.. Let 
A be subspace of L 2 (M") of functions with support in D. Then A ± is 
subspace of L 2 (M. n ) functions with support in D c . Decompose ip(x,tx) as 

with ifii G A and ip2 £ A- 1 . The probability pi of finding particle in D is 



If the result of the measurement at time t\ is that the particle lies in D, 
the variational principle replaces {x(t)} with {£(£)} and, consequently, 
replaces ip(x,ti) = tjJi( x ) + ^2(30) with ip(x,ti) where 



Postulating the variational principle of the previous section (rather than 
the invariance of the phase at ti), we have therefore obtained the so-called 
"collapse of the wavefunction", see e.g. and references therein. The 
collapse is instantaneous, precisely as in the orthodox theory. It occurs 
"when the result of the measurement enters the consciousness of the ob- 
server" |17j. We mention here that, outside of stochastic mechanics, there 
exist alternative stochastic descriptions of (non instantaneous) quantum 
state reduction such as those starting from a stochastic Schrodinger equa- 
tion, see e.g. [3B] and references therein. 

7 Conclusion and outlook 

We shall show elsewhere jlH] that the variational principle of section El 
may be replaced by two stochastic differential games with real velocities 



px(x) 



Xp(x)p(x,t 1 ) 
Id P( x ', t^dx' 




^(x,ti) = p{x)l exp[-S(x,t 1 )} 



V*iQe) 
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with an appealing classical interpretation. We shall also show that, using 
Nelson's observation in (3B1E3 and this variational principle, it is possible 
to obtain a completely satisfactory classical probabilistic description of 
the two-slit experiment. 

If the variational mechanism described here can be extended to the 
case where both the initial and final quantum states are varied, it would 
provide a general approach to the steering problem for quantum systems 
(extending [50]) that has important applications in quantum computa- 
tion [ST], control of molecular dynamics [52] and many other fields. 
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